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I. INTRODUCTION 



Energy levels of classically chaotic systems exhibit correlation only slowly subsiding with the energy offset [T]. 
Around 1980 it became clear, after extensive numerical experiments, that with as few as two degrees of freedom, 
spectral correlations in the highly excited energy domain have universal properties and obey the same laws as the 
eigenvalues of the Gaussian random matrix ensembles of the appropriate symmetry class [2H5] ■ This assertion, known 
as the Bohigas-Giannoni-Schmitt (BGS) conjecture took a surprisingly long time to be proven. In most cases the 
tool used was the Gutzwiller formula giving the spectral density p [E) in a chaotic system as sum over the classical 
periodic orbits, each orbit creating a contribution ~ ^^s^/h -^^j^gj-g jg ^^j^g action of the orbit 7. Its substitution into 
the spectral correlation function and its Fourier transform K (r) , the spectral form factor, leads to double sums over 

orbit pairs with summands proportional to e^{^'^~^'''^^^ . Significant contributions can be expected only from pairs 
with the action difference not large compared with h. 

The first success in the proof of BGS was connected with the diagonal approximation [6 which takes into account 
only pairs with = Sy; it explained the fact that at small times K (t) w 2r (time reversal allowed, orthogonal 
universality class) or t (time reversal forbidden, unitary class). Fifteen years later came the realization that a long 
periodic orbit with a small-angle self-crossing dividing it into two pieces, has a " partner" orbit with the crossing 
avoided, but otherwise almost unchanged, up to the sense of traversal of one of the pieces [7]. Contributions of such 
" Sieber-Richter pairs" sum up the next-to-leading term — 2t^ in the form factor. Summation over pairs in which the 
partner consists of pieces of the original orbit reconnected in all thinkable ways, restores the small-time form factor 
to all orders 

The form factor experiences a break-up of analyticity at t — 1, i.e., at the Heisenberg time Th = 27rfip, which 
reflects the existence of an oscillatory component of the correlation function with the period of the mean level spacing 
A = 1 /p. That component is overlooked in the straightforward semiclassical approach providing correlation functions 
as asymptotic power series in 1/e; the reason is that semiclassical sums need for their convergence a non- vanishing 
positive imaginary part of the energy parameters, however e*^ • -\- • -I- ... if Ime > S > 0. Early estimates 
of the oscillatory components are contained in |10l lllj . Systematic approach is based on the formalism of generating 
functions, i.e., averaged ratios of the spectral determinants. In this approach partnership of more than two classical 
orbits is taken into account, and use is made of the semiclassical approximation of the spectral determinant known as 
the Riemann-Siegel look-alike p]^HT5] . As a result complete agreement of the semiclassical correlation functions with 
RMT for spinless systems was demonstrated [El [17] . 

Systems with half-integer spin belong to the symplectic universality class whose RMT counterpart is the Gaussian 
symplectic ensemble. The spin coupling to chaotic translational motion leads to randomization of the spin evolution 
[18j . The ergodicity of that evolution is instrumental for the evaluation of the relevant periodic orbit expansions such 
as the diagonal sum for the form factor jl9j and the contribution of the Sieber-Richter pairs which changes its sign in 
the presence of a half-integer spin |20| . The full expansion of the form factor of systems with symplectic symmetry 
for times smaller than Th was obtained for the quantum graphs in [21' and for general dynamical systems in [5 . 

Here we close the gap in the proof of BGS for the systems with half-integer spin by showing the equivalence of 
their semiclassical correlation function with the RMT predictions including the oscillatory terms; the corresponding 
form factors coincide for all times. Analytical properties of the correlation function are used to recover the oscillatory 
term with smaller frequency responsible for the well-known logarithmic singularity of the form factor. The derivation 
employs the simple duality discovered between the semiclassical 4-determinant generating functions of chaotic systems 
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with symplectic and orthogonal symmetry. We talk mostly about the spin 1/2 case although the results remain true 
for other half-integer spins. 



II. COMPLEX CORRELATOR, FORM FACTOR AND GENERATING FUNCTION 

Our main object of study will be the complex two-point spectral correlation function (complex correlator, for short). 
This is an analytic function of the complex dimensionless variable e which is the double spectral sum, 

/ 1 1 \ 1 

where A stands for the mean level spacing; (. . .) denotes averaging over an interval of the reference energy E, classically 
small but large compared with A. The complex correlator is defined in the half plane Ime > where it is analytic 
and tends to zero when |e| oo, and can be continued to the lower half plane where ([T]) would no longer be true. The 
real part of C (e) at the positive real axis coincides with the real level-level correlation function [1] while its Fourier 
transform is the spectral form factor K(t)] the connection between the two functions is given by 

if(r) =^(C) = ^/_°l+;°„C(e)e-2-/-de, r > 0; (2) 
C{e) =T-\K)=2j^e''^^^/^"K(T)dT, Ime > 0. (3) 

Here th stands for 1 for the orthogonal universality class and 2 for the symplectic class; such choice is equivalent to 
the replacement 2e — e in the symplectic case motivated by the Kramers degeneracy 

The semiclassical evaluation of the complex correlator is based on the generating function defined as the averaged 
ratio of four spectral determinants, 

^ / det(g-£;-^)det(H-i;-gj^) \ 
W{H-E-^)det{H-E-^^) ■ 



where (e) = (e^, es.eci £n)-The complex correlator can be obtained from Z as 
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Cie)= \iui 2 „ Z 

SB ,£z?— f — e 



(5) 



III. SEMICLASSICAL GENERATING FUNCTION FOR ORTHOGONAL SYMMETRY CLASS 



Here briefly we recapitulate the results for systems with orthogonal symmetry. The semiclassical representation of 
the generating function follows from the chain of relations [121 E] , 



det [H - E) ^ exp 



- / dE' Tr {H ~ E') ^ 



exp 



,iS-,{E)/h 



the last step is the Gutzwiller expansion of the exponent into a sum over periodic orbits 7 with the actions S-^ and 
stability coefficients Expanding all four exponentials we get a sum over quadruplets of " pseudo-orbits" A, B, C, Z?, 



(1) (0\ - 



Z^^^ is) 




FaFcF*sF*o (-1) 



(6) 



^^iAS/h^i{TAeA+Tcec-TBeB-TDeD)/TH 



a pseudo-orbit, say A, is a set of va periodic orbits whose actions and periods sum up to Sa,Ta and whose product 
of stability coefficients is Fa- The difference of actions AS* = Sa + Sc — Sb — Sd must be small compared with h 
for the quadruplets making meaningful contributions. Note the sign factor (—1)'''^"'"'^° depending on the number of 
orbits in the pseudo-orbits associated with the numerator of the generating function. 
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The leading contribution is created by the diagonal quadruplets in which the pseudo-orbit pair {B, D) contains the 
same periodic orbits as (A, C) and consequently AS* = . The diagonal contributions sum up to 

7^ _ ,MeA-eB-ec+eo) (^^ ^ (^C - Spf 

-^diag — e , ^2 / x2 y I 

{£A-eB) {£c-£d) 

and can be factored out like Z^^^ — -^diag (1 + ^off); the off-diagonal part stands for a sum similar to (joj) but with 
all orbits of {A^C) different from those of {B,D). Contributions with a small action mismatch now come from the 
quadruplets in which the periodic orbits of {B,D) are " partners" of those of {A,C) , i.e., consist of practically the 
same but differently connected pieces. Reconnections occur in " encounters" which are places of close approach of 
I > 2 almost parallel stretches of the same or different orbits; the possibility of such a reconnection is a fundamental 
property of chaotic motion [22 . Using ergodicity summation over all pseudo-orbit quadruplets could be reduced to 
summation over their topological families (" structures"). In the end the expansion Zos = ^n, was obtained 

[TUIlT]. Here Zn ^ accumulates contributions of quadruplets with L — V = n; V is the number of encounters 
containing L — stretches. The explicit expression of Z„ for the orthogonal class is (" O" =orthogonal) , 

7 (0\ (£A - £c) {sb - ep) 
(ea - £d) [Sb - ec) 



1)!2" / 1 



n 



The true high-energy asymptotics of the generating function is not exhausted by Z'^^\ Recovery of the missing 
component was achieved on the basis of the so called " Riemann-Siegel look-alike" representation [T^ [T31 [H] of 
the spectral determinant, due to the basic quantum mechanical symmetry properties. The additional summand is 
obtained from Z^^^ (e) by the "Weyl transposition" w (e) = (e^, sb.sd, £c) of its arguments, , 

Zo (e) - ^(''^ = (e) + (e) , (9) 
Z(2)(e) = z(i)(ii;(e)). (10) 

As shown in [17,, Z^^^^^ can be summed up and then yields the RMT generating function of the orthogonal ensemble 
Zgoe- 



IV. SYSTEMS WITH SPIN 1/2. DUALITY WITH SPINLESS CASE 

The spin-1/2 evolution must be treated quantum mechanically while the orbital motion still allows semiclassical 
description. We assume that the spin is driven by the interaction with the translational motion while neglecting the 
back reaction of the spin |18| . The van Vleck propagator of the two-component wave function then falls into a product 
of the semiclassical propagator of a spinless particle, and a 2 x 2 matrix of the spin evolution. The Gutzwiller formula 
follows from the van Vleck propagator after going to the energy representation and taking a trace; consequently 
contribution of a periodic orbit 7 in the Gutzwiller expansions is now to be multiplied by Tr Uj where is an SU2 
matrix describing the change of the spin state after a single traversal of 7 [15] . 

The contribution of a quadruplet {AC){BD) in ([6]) is a product of Gutzwiller amplitudes of all the orbits constituting 
the quadruplet. Therefore, in the presence of spin it has to be multiplied by the product of traces of the spin evolution 
matrices for all its orbits, 2ac,sd = n7e(AC) '^''^ ^jYl-y'eiBD) TrC/^'. We recall that the orbits in (AC) and {BD) are 
constructed of the same L pieces. Therefore due to the group property of the propagator, U^, Uj> can be replaced by 
products of the matrices Vi describing the spin evolution after traversal of the i— th piece. Each 2?^, i = 1, . . . , L, 
occurs once in (AC), once in (BD); if the sense of traversal of the i— th piece is reversed in the partner, the second 
entry would be 

The factor 2ac,sd has to be averaged over an interval of the reference energy E. Interaction with the chaotic 
translational motion makes the spin evolution ergodic [18 . Assuming that Vi associated with different non-overlapping 
pieces are independent quasi-random SU2 matrices, we can thus replace averaging by the integration over all Vi over 
the group SU2- The result for the parthership of just two orbits (the only one of interest at the sub-Heisenberg times) 
is well-known p0l[2T]. 
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where V, L are the number of encounters and encounter stretches in the orbit pair. We generahze it to the pseudo-orbit 
quadruplets containing an arbitrary number of orbits. The most important new element is however that the orbit 
number can change after reconnection in the encounters; see the elementary example Fig. [T] w here two pieces of the 
figure-8 orbit become two separate orbits after reconnection. Calculations in Appendix |VIII B| show that 



\^AC,BD) — 



(-1) 



L — V + VB+I^D — l^A — l^C 



(11) 



where V,L are now the number of encounters and the encounter stretches in the quadruplet (A,C) , (B, D). The 



reflects the difference in the number of orbits after the reconnection {A, C) 



additional sign factor (- 
{B,D). 

The factor (11) leads to important consequences. Namely, inserting it into the semiclassical expansion ^ of the 
generating function and combining with (—1)'^^+'^" we obtain (^_\'^^^^+'^b+vd ^ replacement A,C ^ B,D in 
the exponent shows that the roles of the numerator and denominator in the generating function are reversed. We 
recall that a quadruplet with L — V — n contributes to Z^^^ in the order e^"; our rescaling 2e — )■ e in the symplectic 
case absorbs 2^~^. Consequently all expansion terms of the symplectic off-diagonal sum Zoff are obtained from their 
orthogonal counterparts (jsl) by the interchange A^ C, B ^ D and the sign change of all arguments, 



-£D,—£A, — £b) ■ 



The same substitution connects the full periodic orbit expansions Z^^^ = ^diag (^1 + Zn^^ , 



.(1) 



.(1) 



-£c, —£d, —£a, 



■ £b) 



(12) 



(13) 



We get an important and remarkably simple relation ( 13 ) between the semiclassical generating functions of systems 
with or without half-integer spin. (Turning from semiclassics to RMT we note that numerous identities of that 
kind between the GOE- and CSE-associated functions are well-known under the name of duality relations.) In 



Appendix VIII C we check that a duality relation analogous to ( 13 ) does exist between the 4-determinant generating 



functions of GOE and GSE. A seeming contradiction arises: whereas the high-energy asymptotic expansion of Zqq-e 



is identical with its semiclassical counterpart Zq^\ the analogous expansion of .^gse is not: It differs from Z^^^' (e) 



(12) 



Zi'' is) 



Zg^-* {w (i)) by an additional elementary summand proportional to 



see (32). The reason is 



purely mathematical, and the missing component of the semiclassical generating function can be recovered by Borel 
summation. We prefer to demonstrate the method on the less cumbersome example of the symplectic complex 
correlator, see the next Section. 



The averaged spin factor for spins different from 1/2 is given below in (30). For half- integer spins the result differs 
by the replacement of 2 in the denominator of ( 11 ) by 25" -I- 1, however this is compensated by the changed mean 



level spacing and Heisenberg time [S]; our equation (13) remains in force. For integer spins the sign of the averaged 
spin factor is always positive while 2S -|- 1 in the denominator is compensated in the way just described, therefore the 
generating function is the same as without spin. 
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V. BOREL SUMMATION; MISSING OSCILLATORY COMPONENT AS "STOKES'S SATELLITE 



Applying d'^SASB to Z^^'f (e) + Z^"^^ (e) and going to the limit in (jsj), we obtain the asymptotic expansion of the 
complex correlator for both symmetry classes, 

c(-)-Ef:^+^"^E^ (14) 

n>2 n>4 

where the non-oscillatory and oscillatory parts are generated by Z^-^^^ (e) and Z^'^^ (i) respectively. The coefficients are 
easily calculated from (|8]) and the duality relation; for both symmetries a2 — —1, while for n > 2, 

(n-3)!(n-l) _(n-3)!(n-3) 
an,o- , Vo- , 

The factorial growth of the coefficients signals that the asymptotic series diverge for all e. 

Suppose we want to restore the analytic functions behind these series by means of the Borel method [57 . The first 
stage would be the term-by-term Fourier transform ^ employing 

^ ^ 2i"(n- 1)! ^ ^' ^ ^ 2i"(n- 1)! ^ 

the resulting series in r and (r — 2) converge to analytic functions. On the second stage the inverse Fourier transform 
produces a closed expression for the complex correlator. This is easily done in the orthogonal case and leads to the 



form factor Kgoe (t) and then to the exact Cqoe, see Appendix VIII A Eq. (22) 



The symplectic case is more interesting. On the first stage we obtain the form factor as 

K = e{T) [t/2 - (r/4) In (1 - t)] 

+ e{T-2)[l- t/2 + (t/4) In (r - 1)] ; 

The first summand has a branch cut (l,-|-c»); to proceed with Borel we need to continue ln(l — r) to all r > 1. 
There are three obvious choices: use the logarithm values In |1 — r| ± itt at the lower or upper lip of the cut, or their 
average ln|l — r|. Only the last option is admissible since the form factor must be real, in view of reality of the 
energy eigenvalues [1]; the factor at (r) will then be if < = r/2 — (t/4) In |1 — t| . The Fourier transform of K with 



this choice produces the exact GSE correlator ( 23 1 whose asymptotic expansion contains an additional oscillatory 
summand, 



n=2 n=2 



incidently, the latter would be generated by the term cx g^(^A-eB)/2 Zqse (32) after apphcation of 

It may seem strange that the back-and-forth Fourier transform recovered, free of charge, the missing oscillatory 
contribution to the correlator. In fact, restoration of an oscillatory term, given an asymptotic power series, is a 
legitimate mathematical tool described in detail in the book |24j; see Berry and coauthors ,25. .26] for further important 
developments. The key idea is that the manner in which coefficients of an asymptotic series tend to infinity contains 
information about the exponentially small terms disregarded in the classical Poincare approach. Such terms become 
oscillatory and of crucial importance when the asymptotics is continued to the anti-Stokes lines. 

Here is the barest minimum of detail on the method. Consider a diverging asymptotic expansion of an analytic 
hmction g (z) , z = x + iy, 

oo 

and assume that in the limit of large n its coefficients tend to 

c„^(n-/3)! (16) 



(v-i) ine-^/zf^ ^ 



AS 



FIG. 2: The exponential "satellite" changes by a jump at the Stokes line (red) and becomes oscillatory at the anti-Stokes (AS) 
lines 



Then 



a) The real positive semi-axis is the Stokes line at which the power expansion (15 1 has all its terms positive and 



therefore maximally dominant with respect to an exponentially small additional component gsD (z); the value of 
9SD (z) changes almost by a jump when the positive semi-axis is crossed; 
b) The subdominant component behaves like 



9SD {z) oc 

c) The imaginary semiaxes are the anti-Stokes lines where gsD (z) becomes oscillatory and comparable to the power 
expansion, see Fig. [2j 

d) Under certain assumptions about the properties of g{z) we have, 



aiz) 



R{z) 



n*(z,q) 

E 



n=l 



7/3-1 



R{z), 



i erf 



2/ 

2x 



-V{z,q) 



(17) 



(18) 



The upper limit of the sum is n* (z, q) ~ Int \\z\ -V q\ where q is of the order unity and otherwise arbitrary. The error 



function erf (cr) = ^ e *^dt in (17) is close to 1 for \z\ large and argz > S > 0, and to —1 for argz < —5 <0; on 



the real axis it is zero. The almost jump-like change of the subdominant component when the real axis is crossed, is 
the Stokes phenomenon. The small real correction 



T]{z,q) 



Fract{|z| -l-g} -f /3 



'-3 



2r 
6x 



is significantly non-zero only close to the a;— axis; its dependence on q compensates that of the sum in (17) . The 



constant v must be deduced from additional information on the function g (z). In particular, if g (z) is real on the real 
axis we must choose ^' = 0; the asymptotics of g (z) contains then oscillatory components at both anti-Stokes lines. 

Let us apply the method to the non-oscillatory part of the symplectic correlator. The terms of its expansion 
a„_g/e" are all positive at the positive imaginary axis of e which is the Stokes line where the power series is maximally 
dominant; the anti-Stokes lines are the real semiaxes of e. The results above are applicable with e = iz, /? = 2 
. According to its definition ([T]), the complex correlator must be real for positive imaginary e. Therefore we must 
choose in (17 1 = such that the oscillatory " Stokes satellite" of the power series must be present on both real 
semiaxes, 



R{e) 
R{s) 



'2i' 
2e 



■ -l-oo. 



e — > — oo. 



(19) 
(20) 
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This is indeed the leading term in the oscillatory part of the symplectic correlator recovered by the Borel method. It is 
subdominant in the upper half-plane away from the real axis and experiences the erf —like approximate discontinuity 
at the positive imaginary axis. 

It is instructive to investigate what happens if we choose to continue log (1 — r) as log (r — 1) ± itt for r > 1 in the 
form factor at the first stage of the Borel summation. The functions obtained by the inverse Fourier transform would 
then differ from the correct complex correlator by the additional terms 

/•OO 



They would cancel (19) at one of the real semiaxes and double its amplitude at the other one, i.e., exactly what we 
would get if we chose = ±1 in ( |17[ ). Therefore the alternative choices of v are equivalent to different continuation 
of the result of the first Borel stage beyond the branch point in the t domain. 

Finally let us convince ourselves that additional oscillatory components do not arise in the orthogonal case. The 
e^^— expansion with the coefficients a„_o has its terms all positive on the negative imaginary semi-axis, i.e. in the 
non-physical half-plane of e where ([T]) is inapplicable and the correlator need not be real. The parameter /3 = 2 is the 
same as in the symplectic case but we must now set e = —iz and the Stokes satellite now behaves like ~ e~*^/e. Let 
us change the phase of e from —1^/2 via to positive values; if the satellite were present with a non-zero amplitude 
it would become exponentially large in the physical region Ime > 0. This is forbidden, and we must choose thus in 



(17) V — —1 which corresponds to absence of the exponential term in the sector ~it/2 + 5 < arge < tt, in particular 



at the real positive semi-axis. 



VI. CONCLUSION 



We studied the generating function, complex correlator and form factor of systems with spin 1/2. Expanding the 
generating function into a sum over periodic orbit quadruplets we showed that in the presence of spin, the terms 
of the expansion acquire an additional sign factor whose effect is to interchange the role of the numerator and the 
denominator of the generating function. As a result, the generating functions of the orthogonal and symplectic class 
turn out to be connected by a simple substitution of their arguments. 

The periodic orbit expansion of the generating function supplemented by the Riemann-Siegel look-alike formula for 
the spectral determinants yields the complex correlator as combination of two asymptotic series in the second one 
multiplied by e*^*^; they are responsible for the form factor at small times and at times larger than the Heisenberg time. 
We demonstrate how the Borel summation reveals in the symplectic case one more oscillatory term cx e'^ associated 
with the logarithmic singularity of the form factor; the origin of that term is clarified by the Dingle-Berry method 
of smart summation of the asymptotic series as a display of the Stokes phenomenon. With the missing oscillatory 
term restored, complete equivalence of the correlation functions of chaotic systems with spin 1/2 and the Gaussian 
symplectic ensemble of RMT is reached. 

The inherent ambiguity in restoration of a function from its asymptotic series is solved on the ground of reality of 
the energy eigenvalues. The same reason is at the heart of the Riemann-Siegel look-alike, such that existence of both 
oscillatory components of the complex correlator in the symplectic case can be traced to unitarity of the quantum 
mechanical evolution. 

There are several possible further developments of the theory. An obvious generalization would be to consider 
parametric correlation in systems with spin 1/2 at times comparable with the Heisenberg time 27J. Away from the 
deep semiclassical limit, system-specific deviations from the universal behavior in systems with half-integer spin can 
be of physical interest. 
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VIII. APPENDIX 
A. RMT complex correlator of the orthogonal and symplectic case 



The complex correlators are conveniently expressed in terms of the functions, 

/■oo ±i{t-z) 

f±{z) = j — dt- 

The integral representations are applicable whenever the integral converges; for all z we have 

/± {z) = e^"- - Ci (z) T * Si (z)] . (21) 

where Ci, Si are the integral sine and cosine. These functions are analytic in the plane of z with the cut along the 
negative real axis. 

The complex correlator of the orthogonal case can be represented in terms of /± as, 

(r) (It — Cqqe + Cqq^, 



Cgoe (e) 


/"OO 

= 2 / 

Jo 




^GOE 


1 

~ "2^ 


+ 2^( 


*-"GOE 


~ 





(22) 



while in the symplectic case 



nOO 

Cgse (e) = 2 / e-^i^GSE (r) dr - cg^^ + 4se + 4'se> 
Jo 



.(1) 1 1 



4'sE = -S[H^+l)/+(^)-l]- 



4'sE = ^(^ + 0e-. (23) 

The components C*-^' in both cases have non-oscillatory asymptotic expansion in powers of at the real positive 
axis while C^^-* oscillates like e*^"^. That follows from the asymptotic representations of f± (z), 

f±{z)^-Y,-—^^<j^{z), Izl^oo, (24) 
fc=o (±*^) 
-7r + (5 < argz < TT-O, (/+); 
— TT + < argz < TT — (5. (/-)• 

The Stokes line where the power expansion is dominant is z = —it, t > 0, for (z) and z — it, t > 0, for /_ (z); the 
subdominant satellite, — 27rie~*^ for (7+ (z) and 2-Kie^^ for cr_ (z), exists in the quadrant to the left of the respective 
Stokes line, see Fig. [3j 

The form factors obtained from the complex correlators by the transformation ([2]) are 



TT 



GOE: 



• GSE: 



Js:(r) = 2r-Tln(l + 2T), r < 1; 
i^(r) = 2-rln^^, r > 1. 



K{T)=^--^-\n\l-T\, r<2; 



K{t) = 1, r>2. 
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FIG. 3: Asymptotics of the functions /+ (e) (a) and /_ (e) (b): exponential term is present left of the Stokes line (red) 



B. Spin factor for quadruplets 

Contribution of the quadruplet {AC){BD) in the periodic orbit expansion of the generating function Q contains in 
the symplectic case the spin factor E.ac,bd = Y[-ye(AC) ^7n7'e(-B-D) ^i' "^hich is product of traces of the periodic 
orbits composing the quadruplet. Non- vanishing contributions to the sum are created by the quadruplets such that 
the orbits of {BD) are constructed from pieces of the orbits in [AC) connected in different order and possibly traversed 
with a different sense. We assume ergodicity of the spin motion and independence of the spin evolution along different 
orbit pieces. Averaging is then done step by step by integration over the evolution matrices associated with the 
orbit pieces; the integration domain is the group SU2- Possible outcomes of a single integration are summed up by 
the relations [3T], 

j dVTv {ADBV) = - ^ Tr {AB-^) , (25) 

j dV Tr {AVBV- ^ ) = ^ Tr Tr (26) 

j dV Tr {AV) Tr {BV) = ^ Tr {AB-^) . (27) 

Here A,B are any fixed SU2 matrices. Applying these rules, e. g., to the structure with = 1,L = 2 shown in Fig. 
[1] we obtain the result, 



{^ac.bd) - j dVidV2 Tr (I?i2?2) Tr (2?i) Tr (P2) = \, 



which is opposite in sign compared with the Sieber-Richter pair [20] . 

We shall find the average of '^ac,bd for all structures extending the inductive method of Bolte and Harrison from 
pairs of orbits to the pseudo-orbits quadruplets. We remind the essence of the method. Consider an orbit pair 7, 7' 
differing in V encounters with L encounter stretches and assume that averaging of the spin factor produces the factor 
C-yy = (-l)^"^/2^-^. Introduce one more orbit 7" differing from 7' by an additional 2— encounter such that the 
pair 7, 7" contains V" — V+\ active encounters with L" = L + 2 stretches. It is then shown using the recurrence 

relations that Cj^^' — —Cjj'/2 = (—1)^ ~^ /2^ ~^ . Starting from a pair without active encounters L = V = 
when the formula is true, and adding 2— encounters one by one we obtain that the result is true for an arbitrary 
number of 2— encounters. Finally, reconnection in any Z— encounter with I > 2 can be reduced to I — 1 successive 
reconnections in 2— encounters in / — I steps, with the factor (—1)^"^ /2^~^ correct on all steps. Indeed, reconnection 



10 



in an /—encounter can be described by a permutation of I elements, however, any permutation can be represented as 
a chain of transpositions of just two elements. 

Let us apply this method to the pseudo-orbit quadruplets. Let F = {A, C) be the initial pseudo-orbit pair, and 
F' — {B, D) its partner differing from F in F active encounters with L stretches. Let F" be a pseudo-orbit pair 
differing from F' by reconnection in a single additional 2— encounter whose stretches belong to some orbit 7 in F and 
to 7' in F' where 7, 7' may differ in an arbitrary number of other encounters. The only new case to be considered is 
that of a 2— encounter with almost parallel stretches such that 7' breaks up after reconnection into two orbits 7p and 
7^' belonging to F", see the example in Fig.fTj 

The averaged spin factors for the quadruplets FF' and FF"can be written as multiple integrals over SU2 with the 
Haar measure. 



C 



Crv" = d{ 



d{l)d{i), 

d{i)d{y;) (7;') 



Here d{"f) ,d (7') denote SU2 integration over matrices associated with pieces of 7 and 7' while d (7^' ) d (7^') indicate 
integration over matrices associated with pieces of the disconnected orbit pair in F". Integration d{. . .) is over the 
remaining variables, same for FF' and FF". In the way of induction, let us assume that 



Cpp' — 



(-1) 



L — V+ur — ^T 



(28) 



where vr — i^A-'r vc ^^nd z/p' =1^3 + 1^0 is the number of periodic orbits in F and F', respectively, and prove that the 
analogous formula will be true for Cpp" . 

Consider Fig. [4]where the orbits 7, 7' and the orbit pair 7" are depicted. In 7, 7' we see a parallel crossing; however 
it is not switched between 7,7', i.e., it is inactive and not counted \n.V = Vpr'- On the other hand, it is activated in 
the pair 77" such that 

Vrp"=y + 1, Lpp"=L + 2, VT"=VT' + l- 

In Fig. [4]the orbit pieces adjacent to the encounter and incorporating parts of its stretches are denoted a, 5, c, d; they 
are assumed to coincide in 7,7', i.e., do not contain any additional active encounters. The associated matrices T^a 
etc. involved in the SU2 integration will be denoted by the same letters for compactness, a = Va etc. (we hope that 
d as the integration variable will not be mixed with d as the differential!) The two links attached to the selected 
2— encounter are denoted Li,L2 hi 7 and ^1,^2 i'^ l' ■ FFnlike a,b,c,d, the links Li and L[ need not coincide; indeed, 
Li can contain any amount of encounters active in FF' such that L^ can contain pieces of all orbits of the original 
quasi-orbit pair F other than 7. 

The SU2 integrals can be written, with dLi being a shorthand for dV]^-^ , etc, 

Cpp' = J d {. . .) dLidL2dL\dL'2 j da db dc dd Tr {LiabL2cd) Tr {L'iabL'2cd) ; 

d{...) dLidL2dL[dL'2 j da db dc dd Tr {LiabL2cd) Tr {L[ad) Tr (X'scb) . 

ah and y = cd as the new integration 



TV- 



Let us transform the integrals over a, &, c, d. In Cpp' we can take the matrices x 
variables writing 



dadbdcdd Tr {LiabL2cd) Tr {L'iabL'2cd) 
J dx dy Tr {LixL2y)Tr {L\xL'2y) 



(29) 



dy Tr 



L2yLALA y-'\L 



Tr 













Li [l 




Tr 


L2 (l 


■)1 



The cyclic invariance of the trace and the relations ( 27 ) and ( 26 ) were used. 



11 






FIG. 4: A 2-encounter belongs to the orbits 7 in the pseudo-orbit F and 7' in the partner pseudo-orbit F'. Reconnection in 
that 2-encounter divides 7' into 7" ,72 creating a new pseudo-orbit F" 



Now let us carry out similar transformations of the integral in Crr" introducing consecutively new integration 
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variables x = ad, y = ch, z = cd. 



da dh dc dd Tr (LiabL2cd) Tr (L'^ad) Tr {L'2cb) 
j dxdhdcdd Tr (^Lixd~^bL2cd^ Tr {L'lx) Tr [L^ch) 
d-^bL2cdLi (L^y^] Tr(L^c6) 

yL2cdLi {L[)-'] Tr(4j/) 



_ 1 

^ 2 

_ 1 

~ 2 

_ 1 

" 4 

_ 1 

~ 4 

= iTr 



dfo dc dd Tr 
dc dd Tr 



d-^c-i^ 



dcdd Tr 



L2cdLi{L[) ^d-^c-yL2) ^ 



- / dz Tr 



Tr 



Comparing with ( 29 ) we sec that 



Crr" — 2^^"^' 



which agrees with ( 28 1 with L ^ L" = L + 2,V ^ V" ^ V + l,vr' vp" = i^r' + ^ ; therefore if ( 28 ) is true for the 
quadruplet TV it will also hold for rr". 

Evidently, the reversed process when two orbits of F' merge into a single orbit in F" after reconnection in a 
2— encounter, also agrees with (28 1; reconnection in an ^— encounter with Z > 2 is reducible to a sequence of reconnec- 



tions in 2— encounters. By repeated activation of the encounters resulting in joining and disjoining the orbits, we can 
construct any pseudo-orbit quadruplet out of an orbit pair for which (28) is known to be correct; hence by induction 
it is true for all quadruplets. 



For the spin different from 1/2, the Bolte-Harrison recurrence relations (25), (26 1, (27) differ by the replacement of 



2 in the denominator by 25'+ 1; for integer spins the sign in (25) is plus. Otherwise the reasoning remains unchanged 
with the result ( 28 1 replaced by 



rr' 



(-1) 



(25 + 1 



\L-V 



-, half- integer S; 



Crr' 



(25+1) 



L-V 



, integer S. 



(30) 



C. Generating functions of RMT. GOE-GSE duality 

The generating function of GOE found by Zirnbauer |28j has a Weyl symmetric form, 

Zgoe (e) = Fgoe (i) + -Fgoe (w (e)) 
where i = (sa^b^c^d)', w interchanges C and D. Denoting ad = ea — £d etc and assuming Ivaab > we can write. 



^GOE (e) = e 



i ^ {ah—cd) 



ad cb 
ab cd 



1 + -ac bd _^ /+ 
^ cd 



with defined above, see (21). 
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The generating function of GSE is given in pS] as, 



ac ad be bd I 

'^GSE {£) — ^GUE (ej H -— — i 7(^-1 

ao erf 4 



Go 



Gi {z)^i e'^iqdq 



Go (x) = 



sm qx 
1 g 



dq = TT 



iq 



-dq 



00 iga: 



-rf(7 



(tt - + [x) - G<,^^ + G^'^. 



Here ^gue (e) is the generating function of the unitary ensemble; in the last line x = crf/2 is assumed real positive. 

The GSE generating function can be transformed to a sum of two components connected by the Weyl substitution 
w. First let us write Zgse (e) = F^^^ + -F"^^^ with 



(ab-cd) ad cb 



ab cd 



1 —(2-iab) , 

i + -ac bd ^ ' j; 

^ ab 



^ arf ch ac bd (2 — i ab) 



and 



i?(2) 



- e' 2 

2a6 erf 



.(ob+cd) ac db 



ab dc 



ab 



1 (2-iab) 

1+ ad bc ^ V + 
2 a6 



/ M \ ■ (ab — cd) 

The part i^'- proportional to e* 2 generates the non-oscillatory component of the symplectic correlator while the 

part proportional to e'"2" generates the correlator term proportional to cx e*^. The part ' oc 2 is responsible 
for the component of the correlator cx e'^^. 

Let us define the function /_|_ {—x) for real positive x, i.e., at the branch cut, as the average of the values of the 
analytic function /+ (z) at the lips of the cut. 



Taking into account that 

/+ i-x + to) = /; (x) , 

f+{-x + iO) - f+{~x - iO) = -27rie", 

we have Gq (x) = — «/+ a^) e^*^ + if^ (x) e'^: we have concealed tt in the first summand. Now considering that the 
substitution w changes the sign of erf we obtain 



^GSE (e) = ^GSE (e) + -Fgse {w (e)) . 



(31) 



with 



Fgse (e) = = e 
The duality relation now holds. 



{ab-ad) ad cb 

2 

■ erf 



bd 



2 — iab 
— 2 
ab 



u - 



cd 



Fgse {sa, £b,£c, £d) — -Fgoe (— ec, —en, — ea, — £s) : 
Re(eA-eB)>0, Im (e^ - e^) = +iO, 
- en > 0. 
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It is instructive to compare the RMT function (31 ) with the semiclassical results. Taking (IS]) into account and using 



the semiclassical duality, we have the following equivalence in the high-energy limit, 



i^GSE (e) ^ (e) + e*(----«)/2 ---- ^^—f^n, (32) 



ad cb ac bd (2i + ab) 



2ab cd q5 



FosE{w{e))^zf^ {i)^Z^P {w{i)). 



We stress that the asymptotics of Fqse {w (e)) is not simply the Weyl-transposed asymptotics of Fqse (e)- The reason 
is the Stokes phenomenon; the Weyl operation changes the sign of cd , and whereas the asymptotics of /+ (cd/2) is 
purely power-like, /+ (— c(i/2) contains an additional exponential summand, in accordance with 



/+ (x) - (T+ (x) , /+ {-x) - iTre" + cr+ {-x) , 

X — > +0O. 



where cr+ is the asymptotic power series defined in ( 24 1 
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